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ABSTRACT 

The nearby Sun-like star HD 114174 exhibits a strong and persistent Doppler acceleration 
indicating the presence of an unseen distant companion. We have acquired high-contrast imag- 
ing observations of this star using NIRC2 at Keck and report the direct detection of the body 
responsible for causing the "trend". HD 114174 B has a projected separation of 692 ±9 mas (18.1 
AU) and is 10.75 ± 0.12 magnitudes (contrast of 5 x 10^^) fainter than its host in the -ftT-band, 
requiring aggressive point-spread function subtraction to identify. Our astrometric time baseline 
of 1.4 years demonstrates physical association through common proper motion. We find that 
the companion has absolute magnitude, Mj — 13.97 ± 0.11, and colors, J — K = 0.12 ± 0.16 
mag. These characteristics are consistent with an wT3 dwarf, initially leading us to believe that 
HD 114174 B was a substellar object. However, a dynamical analysis that combines radial ve- 
locity measurements with available imaging data indicates a minimum mass of 0.260 ± O.OIOMq. 
We conclude that HD 114174 B must be a white dwarf. Assuming a hydrogen-rich composition, 
atmospheric and evolutionary model fits yield an effective temperature Toff — 8160 ± 4000 K, 
surface gravity log 5 = 8.90 ± 0.02, and cooling age of tc ~ 3-4 Gyr, which is consistent with the 
4.7^2 I Gyr host star isochronal age estimate. HD 114174 B is a benchmark object located only 
26.14 ± 0.37 pc from the Sun. It may be studied at a level of detail comparable to Sirius and 
Procyon, and used to understand the link between the mass of white dwarf remnants with that 
of their progenitors. 

Subject headings: keywords: techniques: radial velocities, high angular resolution; astrometry; stars: 
individual (HD 114174), white dwarfs 

1. INTRODUCTION 
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The TRENDS (TaRgetting bENchmark-objects 
with Doppler Spectroscopy) high-contrast imag- 
ing program uses precise radial velocity (RV) 
measurements obtained over a long time base- 
line to identify promising targets for follow-up 
observations using adaptive optics (AO) (Crepp 
et al. 2012b, Crepp et al. 2012c). In ad- 
dition to diffraction-limited imaging, TRENDS 
also uses coronagraphy, point-spread function 
(PSF) subtraction, and «1 hour integration times 
to achie ve deep contrast at near-infrared wave- 
lengths dMarois et aLll2006l: iLafreniere et al.ll2007 : 



Crepp et al.ll2010f) . The primary goal of TRENDS 
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is to image and characterize benchmark M-dwarfs, 
brown dwarfs, and (ultimately) extrasolar planets, 
using a combination of complementary substellar 
companion detection techniques, in order to cal- 
ibrate theoretical evolutionary models and theo- 
retical spectral models of cool dwarf atmospheres. 

Precise Doppler measurements for TRENDS 
are obtained at Keck Observatory using the 
High Reso l ution Echelle Spectrometer (HIRES; 



Vogt et alj (|1994l )) as part of the Califo rnia Planet 



Search program ([Howard et alJ l2010t ). Targets 



with the longest RV time baselines have comple- 
mentary data taken at Lick Observatory using 
the 0.6m Coude Auxiliary Telescope (CAT) and 
Shan e 3m Teles cope with the Hamilton spectrom- 
eter (|Vogtlll987l ). High-contrast imaging follow-up 
observations are conducted at Keck using NIRC2 
(instrument PI: Keith Matthews) with the Keck 



II AO system ( Wizinowich et al. 2000l) . and 



at 



Palomar using Project 1640 with the PALM - 3000 
high-order AO system ( Hinklev et al. 201 if ) . 
CrePD et all l2012bl (hereafter TRENDS I) re- 



ported the discovery of three benchmark high- 
mass ratio binaries, HD 53665, HD 68017, and 
HD 71881, each companion having an age and 
metallicit y estimate inferred from the solar-type 
primary ( Crepp et al. 2012bfl . Arguably the most 
interesting object of the three, HD 68017 B, has 
a RV time baseline of 15 years and shows signif- 
icant orbit curvature. With a projected separa- 
tion of only 13.0 AU, it will be possible to mea- 
sure a precise dynamical mass in the next several 
years. The mo del- dependent mass of HD 68017 B 
from relative photometry is rTimodoi — 168 ± 21Mj 
(0 .16±0.02Mq). 

ICrepp et al.ll2012d (TRENDS II) reported the 
direct detection of an TOmodei = O.52±O.O4M0 ter- 
tiary companion orbiting the 84 day period, single- 
line spectroscopic binary HD 8375 AB. With a 
separation of 0'.'3, HD 8375 C has a brightness 
comparable to residual scattered starlight from the 
inner binary pair prior to PSF subtraction. The 
HD 8375 system represents one of many hierar- 
chical triples in the solar neighborhood that are 
known to exist, through statistical analysis of mul- 
tiplicity fractions for stellar systems wit hin < 10 
pc, b ut so often evade direct detection ([Tokovinin 
\2004) . 

In this paper (TRENDS HI), we report the 
discovery of a faint object orbiting the nearby, 



d = 26.14 ± 0.37 pc (jvan Leeuwenll2007f) . G5 star 
HD 114174 (LHS 2681, HIP 64150, GJ 9429). 
The companion, HD 114174 B, cannot be seen 
without significant on-source integration time 
and PSF subtraction. It has a projected sep- 
aration (0.69" — 18.1 AU) comparable to the 
semimajor axis (a = 18.3 AU) of the brown 
dwarf HR 7672 B, which recently yielded a dy- 
namical mass with a fractional error of only 4% 
(jCrepp et al. I l2012al) . Initially characterized as a 
substellar object based on its brightness and col- 
ors, we find that HD 114174 B must instead be a 
white dwarf, in order to establish a self-consistent 
interpretation with long-term Doppler RV mea- 
surements. 

Few "Sirius-like" systems, i.e., directly imaged 
white dwarfs around solar-type stars, are presently 
known - appro ximately one dozen within 20 pc 
(iHolberdliofih . Our serendipitous discovery of 
a compact object orbiting HD 114174 thus con- 
stitutes a new member of a small collection of 
directly detected benchmark white dwarfs for 
which we can determine physical properties in- 
dependent of spectro-photometric measurements. 
For example, other than Sirius and Procyon, 
HD 114174 B appears to have the shortest or- 
bital period a mongst nearby Sirius-like systems 
(iHolberd l2009f ). Thus, forthcoming measurements 
will ultimately lead to tight dynamical constraints 
on its orbit and mass. More than 16 years of RV 
measurements are already in hand. 

2. OBSERVATIONS 

2.1. High-Resolution Stellar Spectroscopy 

Doppler Measurements 

We obtained RV data of HD 114174 using 
HIRES at Keck. The standard iodine cell ref- 
erencing method was used to calibrate instru- 
ment drifts and measur e precise Doppler shifts 
(jMarcv and Butleiill992r ). Our Doppler observa- 
tions began on January 14, 1997 UT. Several years 
of measurements revealed that the star exhibits a 
(mostly) linear acceleration indicative of a distant 
companion. Fig. 1 shows Doppler observations 
taken over an 16.0 year timeframe. RV values, 
epochs, and measurement uncertainties are listed 
in Table 1. 

Treating the Doppler trend as purely lin- 
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Fig. 1. — Relative radial velocity measurements of HD 114174 from the years 1997 through 2013. (Top) 
A long-term Doppler acceleration of 61.1 ± 0.1 m/s/yr implies the presence of a companion with a long 
period. (Bottom) A linear fit model results in large residuals indicating systemic curvature (change in the 
acceleration). Thus, it will soon be possible to calculate the three-dimensional orbit and dynamical mass of 
HD 114174 B. 



ear and using a Markov Chain Monte Carlo 
fitting procedure, we find an acceleration of 
61.1 ± 0.1 ms~"'^yr~^. This analysis includes a 
differential offset between measurements obtained 
before and after the summer 2004 HIRES detec- 
tor upgrade, which we treat as a free parameter. 
Closer inspection of the data shows that the RV 
slope is slowly changing with time. For instance, 
the bottom panel of Fig. 1 displays RV residu- 
als once the measurements are subtracted from a 
linear model. Continued Doppler monitoring will 
allow us to place tight constraints on the com- 



panion orbit and dynamical mass, as soon as the 
direct imaging astrometry shows curvature. 

Stellar Properties 

Stellar template spectra, taken with the io- 
dine gas cell removed from the optical path, 
were analyzed using the LTE spectral synthesis 
code Spectroscopy Made Ea sy (S ME) described in 
Valenti and Piskunov (Il996l) and lValent i and Fischeii 
(|2005l ) . SME provides an estimate of the stellar ef- 
fective temperature (T^ff), surface gravity (log^), 
metallicity ([Fe/H]), and projected rotational ve- 
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locity (iisini). The estimated physical properties 
of HD 114174 derived from spectral fitting are 
shown in Table 3. 

HD 114174 is listed in the SIMBAD database 
as a G5(IV) subgiant. However, its location with 
respect to th e mean Hipparcos main-sequence 
(IWrightl [20051) and spectroscopic \ogg = 4.51 ± 
0.06 value indicate the star is a main-sequence 
dwarf. Placing HD 114174 on a Hertzsprung- 
Russell diagram, we estimate a mass of 1.05 ± 
O.O5M0 and age of 4.7^2^5 comparing 
spectroscopic parameters to Yonsei-Yale theoret- 
ical isochrone tracks. Despite its close distance, 
HD 114174 would rarely, if ever, make typical 
high-contrast imaging target lists, which nom- 
inally only select stars younger than «1 Gyr 
( Crepp and JohnsonI 2011 ). 



2.2. High-Contrast Imaging 

We acquired first-epoch high-contrast images of 
HD 114174 on 2011 February 22 UT using NIRC2 



and the Keck II AO system (IWizinowich et al 
2000|) . The bright {Ks = 5.20. ISkrutskie et al 



20061 ) star was placed behind the 300 mas diameter 
coronagraphic spot. We used the angular differen- 
tial imaging (ADI) technique to discriminate be- 
tween residual scattered starlight and companions 
through point -sprea d- function (PSF) subtraction 
jMarois et^l2006l ). 

Raw frames were processed by flat-fielding the 
array, replacing bad pixels with interpolated val- 
ues, and high-pass Fourier filtering. The latter 
step reduces the AO halo and thermal background 
(pedestal) while accentuating quasi-static speck- 
les, which aids with image alignment. Frames were 
then multiplied by a gray-scale intensity mask 
(having contiguous transmission values between 
and 1) to prevent the coronagraph from biasing 
image alignment. We used the Locally Optimized 
Combination of Images (LOCI) algorithm to in 



creas e effective contrast levels (jLafreniere et al 
20071 ). The subtracted frames were then rotated 



to the standard East-North orientation using the 
parallactic angle to combine images and stack light 
from any off-axis sources. 

Upon processing the full 2011 February data 
set, we noticed a faint candidate companion lo- 
cated ss 0.7" to the south in the final reduced im- 
age. Follow-up observations were taken on 2012 



HJD-2,450,000 


RV [m s-'-] 


Uncertainty [m s ^] 


463.1438 


-403.39 


1.27 


545.9798 


-390.77 


1.29 


787.1381 


-351.45 


1.89 


787.1439 


-337.14 


1.33 


838.0563 


-333.70 


1.17 


838.1109 


-329.74 


1.10 


839.1258 


-334.99 


1.20 


862.0168 


-319.81 


1.06 


862.9814 


-320.76 


1.30 


955.8422 


-313.57 


1.45 


981.7798 


-311.44 


1.18 


982.7957 


-312.89 


1.09 


1009.8359 


-302.12 


1.13 


1050.7552 


-288.92 


1.92 


1050.7572 


-290.58 


1.29 


1171.1521 


-273.81 


1.35 


1172.1624 


-273.77 


1.22 


1173.1611 


-277.47 


1.31 


1227.0404 


-267.37 


1.28 


1228.0764 


-259.05 


1.17 


1228.9999 


-262.77 


1.37 


1310.8726 


-249.58 


1.33 


1311.8551 


-243.23 


1.27 


1312.8933 


-247.37 


1.42 


1313.9960 


-249.02 


1.40 


1340.8177 


-242.25 


1.39 


1341.8676 


-248.12 


1.49 


1342.8364 


-242.09 


1.21 


1553.1269 


-205.71 


1.55 


1553.1298 


-200.00 


1.22 


1582.0826 


-196.36 


1.34 


1583.0577 


-197.94 


1.41 


1585.0085 


-196.73 


1.20 


1585.0130 


-196.74 


1.19 


1586.0218 


-205.15 


1.34 


1679.8938 


-179.28 


1.42 


1702.9103 


-175.20 


1.28 


1703.8065 


-174.94 


1.35 


1704.8793 


-177.12 


1.49 


1705.8784 


-179.77 


1.35 


1706.8757 


-175.24 


1.43 


1755.7542 


-168.56 


1.11 



Table 1: 

February 02, 2012 May 29, and 2012 July 4 UT 
with the K' and J filters to obtain color infor- 
mation and determine whether the candidate is 
associated with the primary. Fig. 2 shows the 
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HD 114174 J, 



0.5" 

NIRC2, K' 

13.1 AU February 2, 2012 



HJD-2,450,000 


RV [m s-'-] 


Uncertainty [m s ^] 


1884.1599 


-141.95 


1.41 


1898.1674 


-141.82 


1.28 


1899.1728 


-137.89 


1.37 


1900.1678 


-139.61 


1.30 


1901.1771 


-138.81 


1.27 


1972.1057 


-129.95 


1.26 


1973.1304 


-127.04 


1.33 


1982.1034 


-126.24 


1.43 


2030.9064 


-125.05 


1.43 


2101.8222 


-117.06 


1.75 


2334.1177 


-67.85 


1.84 


2712.0408 


-0.36 


1.37 


3071.9391 


54.34 


1.56 


RV Offset 


3480.0014 


130.10 


0.86 


3480.0026 


131.17 


0.62 


5232.1301 


414.79 


0.91 


5973.1039 


529.01 


0.95 


6109.7547 


554.22 


0.92 


6319.0020 


584.50 


1.10 



Table 2: Doppler radial velocity measurements for 
HD 114174. A horizontal line indicates the divi- 
sion between HIRES (detector) pre-upgrade and 
post-upgrade data points, requiring a differential 
offset for analysis. Uncertainties correspond to 
photon-noise alone. 



confirmation image taken one year after the first 
epoch. 

We also obtained a deep ADI sequence of 
HD 114174 with NIRC2 in the L'-band on 2012 
May 07 UT. These observations were recorded un- 
der good conditions (RiO.5" seeing), at low airmass 
(1.03-1.05), and achieved 40.5 degrees of field ro- 
tation, yet resulted in a non-detection, suggesting 
that HD 114174 B may have blue colors. We derive 
a lower-limit to the brightness difference between 
HD 114174 B and its parent star of AL' > 9.84 
mag. The companion was recovered at all other 
epochs using filters at shorter wavelengths. 

3. PHOTOMETRY 

PSF subtraction modifies the measured fiux 
(photometry) and location (astrometry) of faint 
companions th at are initial l y bur i ed in bright 
speck le noise (jCrepp et al. 1 boill : IPuevo et al.l 
20121 ). We injected synthetic companions into 



Fig. 2. — Confirmation image of HD 114174 B 
taken with NIRC2 at Keck. The AK' = 10.73 ± 
0.12 companion cannot be seen in individual raw 
frames with long exposures, but is clearly visible 
after PSF subtraction. 



the NIRC2 ADI data to estimate the amount of 
companion flux lost as a result o f using the LOCI 
algorithm (jSoummer et al.l 120121) . Self-consistent 
introduction of fake companions into each frame 
also allows us to measure the intrinsic scatter in 
photometric and astrometric quantities and thus 
compute their uncertainty. 

We measured the relative brightness between 
HD 114174 B and its host star by comparing the 
fiux from LOCI-processed images to unocculted 
images of the primary star, accounting for the dif- 
ference in integration time. A fast-readout, subar- 
ray mode prevented HD 114174 A from saturating 
the array in unocculted frames. Subsequent reduc- 
tions using fake companions and the same LOCI 
parameters were then performed to determine the 
(throughput) response of the data pipeline. The 
injected companions were placed at approximately 
the same angular separation as HD 114174 B, but 
with azimuthal separati ons sufficient to pr event 
artificial self-subtraction (jOuanz et al.ll2012r) . Un- 
occulted frames of the primary star informed our 
choice of the FWHM for fake companions. 

Measurements using two different filters, K'^^jj^q 
and JmkOi provided relative fiux and color in- 
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HD 114174 Properties 



right ascension 


[J2000] 


13 08 51.02 


declination [J2000] 


+05 12 26.06 


B 




7.47 


V 




6.80 


<^2MASS 




5.613 ±0.026 


Ks 




5.202 ±0.023 


Jmko 




5.584 ±0.027 


Kmko 




5.189 ±0.024 


d [pc] 




26.14 ±0.37 


proper motion 


mas / yr] 


84.72 ±0.59 E 






-670.11 ±0.47 N 


Host Star 


Mass [Mq] 




1.05 ±0.05 


Radius [Rq] 




1.06 


Luminosity [Lq] 
Age [Gyr] 
[Fe/H] 


1.13 

4 7+2-3 
^■'-2.6 

0.07 ±0.03 


log g [cm s 2] 




4.51 ±0.06 


Teff [K] 




5781 ± 44 


Spectral Type 




G5 IV 


V sini [km/s] 




1.8 ±0.5 




Companion 


a 


AJ 




10.48 ±0.11 


AX 




10.75 ±0.12 


AL' 




> 9.84 


J 




16.06 ±0.11 


K 




15.94 ±0.12 


L 




> 14.99 


Mj 




13.98 ±0.11 


Mk 




13.85 ±0.12 


Ml> 




> 12.90 






8160 ± 4000 


log 5 [cm/s^] 




8.90 ±0.02 


tc [Gyr] 

TOdyn [Mq] 
TOinodcl [Mq] 




« 3.4 
> 0.260 ±0.010 
1.15 ±0.01 



Table 3: (Top) Coordinates, apparent magnitude s, distance, and prop er motion of HD 11417 4 from the 



SIMBAD database. Magnitudes are from 2MASS (S krutskie et al.ll2006n . Filter conversions from lCarpente: 



([2001) have been applied to estimate stellar M KO values. The par allax-based distance is from Hipparcos 
measurements using the refined data reduction of van Leeuwen 20071 (Middle) Host sta r physical properties 
are e stimated from SME using HIRES template spectra as well as theoretical isochrones (jValenti and Fischer 
20051 ). (Bottom) Companion magnitude difference, apparent magnitude, absolute magnitude, and compan- 
ion physical properties, such as the estimated cooling age (tc), mass constraint (lower-limit) from dynamics 
(TOdyn), and mass (mmodci) from photometry using white dwarf evolutionary models assuming a pure hydro- 
gen atmosphere. 



formation. Optimized for the local atmospheric 
transmission and thermal background on Mauna 



Kea (|Tokunaga et al.l l2002f ) 

our filter choice for search-mode operation. We 



the K^jj^Q band is 
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measure flux ratios of Aifj^^^Q — 10.73±0.12 and 
1\Jmko = 10.48 ± 0.11 (Table 3). We then con- 
vert the J2MASS and Kg magnitudes of the primary 
star to Jmko and Kmko using the J2MASS - 
color fro m SIMBAD and empirical transforma- 
tions from ICarpenteii (|200l[ ) . 

To place the companion flux in the same fil- 
ter basif0, we downloaded high-resolution spec- 
tra of various cold dwarfs from the IRTF spec- 
tral librar5nas well as the filter profiles from Keck 
Observator5|f|. IRTF library spectra from a G5 
dwarf, HD 165185, served as a fiducial to com- 
pare the colors of the companion to the G5 pri- 
mary, HD 114174 A. HD 165185 and HD 114174 A 
have J — Ks colors consistent to within 0.045 



mag (|Skrutskie et al.ll2006l ). smaller than the 0.12 
mag uncertainty in the measured companion-to- 
star fiux difference. 

Our initial interpretation of the companion's 
properties led us to believe it was a brown dwarf. 
We thus began filter conversion calculations using 
IRTF spectra f r om tw o cold brown dwarfs listed in 
ICushing et all (l2005h . SDSS J125453.90-012247.4 
(T2) and 2MASS J05591915-1404489 (T4.5). We 
first scaled the flux of the different T-dwarf spectra 
by a multiplicative factor to match our I^K'j^j^q 
measurement, then compared the IS.Jmko model 
result to our A.Jmko measurement. The cooler 
T4.5 dwarf spectrum resulted in a significantly 
better fit to our measurements {IS.Jmko within 
la) than the T2 dwarf spectrum which was dis- 
crepant at the 5(7 level. Using the T4.5 spectrum 
reference, we find /S.Kmko — 10.75 ± 0.12 (from 
the measured AK'^jj^q — 10.73 ±0.12), Kmko = 
15.94±0.12, and { J-K)mko = 0.12±0.16 (Table 
3). 

Both the absolute magnitude and colors of 
HD 114174 B are consistent with a cold brown 
dwarf (Fig. 3). However, as shown below (Section 
5.1, Section 5.2), our dynamical analysis precludes 
this interpretation. Since the spectral energy dis- 
tribution of cool white dwarfs and late-type brown 
dwarfs are both often dominated by coUisionally 
induc ed H2 absorption (jBorvsow and Frommhold 
I990I ). and the correction factor between the 



^MKO ^"^"^ I^MKO filters are 
: 2.196 iim. 



^Central wavelength for the _ff',,, 

Ajf/ = 2.124 fim and Xjf 
^ http: / / irtfwob. ifa. hawaii.edu /~spex /IRTF_Spectral_Library / 

http: / / www2 . keck, hawaii . cdu / inst /nirc2old/ filters .html 
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iLeggett et al 



brown dwarfs from 
Black data points represent 
derived from moderate resolution 



(I2OIOI I. 

spectral types 
spectroscopy 
(http://www.gemini.edu/staff/sleggett). Our 
la measurement uncertainties from broadband 
photometry are shown as horizontal lines. 
HD 114174 B appears to be consistent with a 
T3.5 dwarf, but dynamical considerations indicate 
otherwise. Given the comparable colors between 
brown dwarfs and white dwarfs, compact objects 
like HD 114174 B might be considered as a type 
of false-positive for high-contrast imaging surveys 
dedicated to searching for substellar companions. 



^MKO ^^'^ Kmko filters is small, particularly 
for dwarfs with neutral colors, we adopt the 
AKmko — 10.75±0.12 value for subsequent anal- 
ysis. 
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4. ASTROMETRY 

Our astrometric observations consist of four 
epochs taken over an 1.4 year baseline (Table 4). 
The proper-motion of HD 114174 is high, 675 mas 
yr~^, allowing us to easily determine whether the 
companion shares the same space motion as the 
primary. We measured an accurate separation 
an d position angle usin g the technique described 



Crepp et al 



(|2012aD. Following PSF subtrac- 
tion, we fit Gaussian functions to the stellar and 
companion point-spread function to locate their 
centroids in each frame. The results from trials 
using numerous fake injected companions are av- 
eraged and uncertainty in the separation and posi- 
tion angle are initially taken as the measurement 
standard deviation. We then account for uncer- 
taint y in the plate sc ale and instrument orienta- 
tion (IGhez et al.ll20 08^ by propagating these errors 
to the final calculated position. Finally, a 5 mas 
uncertainty is added in quadrature to the rectilin- 
ear position of the star to account for unknown 
distortions introduced by the focal plane occult- 
ing mask (Q. Konopacky 2011, private communi- 
cation). For companions with intensity compara- 
ble to, or fainter than, stellar quasi-static speckles 
using (near) Nyquist-sampled data, it is difficult 
to measure an astrometric position with precision 
much less than a single pixel (Marois et al. 2010t 
Crepp et al.1l20lll ). 

Fig. 4 shows multi-epoch astrometry mea- 
surements plotted against the expected mo- 
tion of a distant background object. We find 
that HD 114174 B is clearly associated with 
HD 114174 A. The position of the companion 
has changed by only 28 mas over an 1.4 year time- 
frame. An unrelated background source placed at 
infinite distance with zero proper-motion would 
have moved relative to the host star by 902 mas 
over the same timeframe. HD 114174 B has a pro- 
jected separation of 18.09 ± 0.34 AU (2012 July 
4 UT) that appears to be decreasing with time. 
Early indications are that the companion sky- 
projected orbital motion is in a counter-clockwise 
direction (north up, east left), but further ob- 
servations are required to determine whether the 
position angle is increasing. 
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Fig. 4. — Astrometric measurements demonstrat- 
ing that HD 114174 B is co-moving with its host 
star. Location of the primary defines the coordi- 
nate system and is denoted by a star. The dashed 
curve shows the path that a distant background 
object would execute from 2011 February 22 UT 
through 2012 July 4 UT, accounting for proper- 
motion and parallactic motion. 



5. COMPANION PROPERTIES 
5.1. Mass Lower-limit from Dynamics 

A lower-limit to the companion mass may be 
determined from the RV accelera tion and di- 
rect imaging angular separation (|Torresl Il999l: 
Liu et albooilHoward et al.ll2010[ ). Although the 
HD 114174 Doppler signal already shows a hint 
of curvature (Fig. 1), we do not compute the full 
range of possible masses because the astrometry 
has yet to show significant orbital motion, which is 
required to place interesting constrain ts on orbital 
elements and the companion mass ([Kalas et al 



20081: ICurrieet al.ll2Qlll:[Crepp et al.ll2012ai bl). In 

stead, we evaluate the local slope of the RV time 
series. We have obtained a recent Doppler RV 
observation (2013 January 26) to facilitate this 
calculation. Using an MCMC simulation to fit 
the instantaneous RV slope, we find a best-fitting 
Doppler acceleration of 57.4 ± 0.4 m s~^yr~'^, a 
value similar to that found in Section 2.1 when 
considering the full RV data set because curvature 
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Date [UT] 


JD-2,450,000 


Filtcr 


p [mas] 


Position Angle [°] 


Projected Separation [AU] 


2011-02-22 


5615.1 


K' 


719.8 ±6.6 


171.8° ±0.5° 


18.82 ±0.32 


2012-02-02 


5960.1 


K' 


701.1 ±5.0 


172.1° ±0.4° 


18.33 ±0.29 


2012-05-29 


6076.8 


K' 


695.8 ±5.8 


170.5° ±0.6° 


18.19 ±0.25 


2012-07-04 


6112.8 


J 


692.1 ±8.7 


172.2° ±0.8° 


18.09 ±0.34 



Table 4: Summary of astrometrie measurements. 



in the RV time-series is still subtle. This measure- 
ment translates to a minimum dynamical mass of 
0.260 ± O.OIOM0, thus precluding the interpreta- 
tion that HD 114174 B is substellar, despite its 
low luminosity and neutral colors. 

5.2. White Dwarf Interpretation 

White dwarfs are a common outcome of the 
evolutionary process for low-mass stars and are 
thus expected to comprise a non-negligible frac- 



tion of stellar co mpanions ( Holberg et al. 1998t 



Breton et al.ll2012V With a sample of more than 
100 targets, it is therefore not surprising that 
the TRENDS program would discover a com- 
pact object orbiting a nearby star. Indeed, it 
seems the only self-consistent interpretation of our 
photometric and dynamical measurements is that 
HD 114174 B is a white dwarf companion. Such 
serendipitous discoveries might be considered as 
direct imaging false-positives when searching for 
exoplanets and brown dwarfs (jZurlo et al.H2013i) . 

Despite the pr evalence of white dwarfs th rough- 
out the galaxy ( Oppenheimer et al] 12001 ). rela- 
tively few b enchmark compact ob jects have been 
discovered ( Muirhead et al. 20131 ) and only sev- 
eral nearby companions with precise parallax mea- 
surements, such as Sir ius B and Prqcyon B, have 
been imaged directly ( Liebert et al. 20051) . Sim- 
ilar to field brown dwarfs, most detailed white 
dwarf studies involve ag e calibration based on stel- 
lar cluster membership (jWilliams et al. Il2004l) . and 
those with reliable ages rarel y have simultane- 
ous dynam ical mass estimates ( Girard et al. 2000l : 
Bondll2009h . HD 114174 thus represents a conve- 
nient system for studying the intricacies of post- 
main-sequence stellar evolution at low mass. 

Without spectra, we cannot determine the 
chemical composition of HD 114174 B. Instead, 
we proceed by considering two different canoni- 
cal cases: pure hydrogen and pure helium white 
dwarf atmospheres. By fitting theoretical models 



to our J and K observations, we can constrain 
the effective temperature, mass, and cooling age 
of HD 114174 B under reasonable assumptions 
about its composition. 

5.3. Synthetic Photometry 8z Fitting 

For pure hydrogen atmospheres and syn- 
the tic spectra, we use t he theoretical models 
of iTremblav et~aLl (|201l[) FI The models in- 
clude improved calculati ons for Stark broaden- 
ing of hydrogen lines ([Tremblav and Bergeron 
20091 ) and non-ideal perturbations, described 
within the occupation prob ability formalism of 
Hummer and Mihalas ( 1988f ). directly inside the 
line profile calculations. The models also in- 
clude Lyman alpha profile calculations from 
Kowalski and SaumonI (|2006l ). Convective en- 
ergy transport is taken into account following 
the revised ML2/a = 0.8 prescripti o n of m ixing- 
length theory from I Tremblav et al. I (I2OIOI) . De- 
tails regardi ng our helium a. t mosph ere models are 
provided in [Bergeron et al. ( 2011 ). Model grids 
cover a range in effective temperature between 
Tcff = 1500 K and 40,000 K in steps varying 
from 500 K to 5000 K. The surface gravity ranges 
from 7.0 < logg < 9.0 by steps of 0.5 dex. Syn- 
thetic colors are obtained us ing the procedure of 
Holberg and Bergeron J 20061) based on Vega fiuxes 
taken from lBohlin and Gillilandl (|2004l ). 

The method used to fit photometric dat a issim- 
ilar to that described in [Bergeron et al.l (|200l[ ). 



which we briefly summarize. We first transform 
the magnitudes in each bandpass into observed av- 
erage fluxes f\,rn using the following equation: 



TO = -2.5 log /a,™ + Cm, 



where 



fx^m — 



(1) 

(2) 



*http: / /www. astro. umontrcal.ca/ ~bergeron/CoolingModels/ 
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Fig. 5. — White dwarf model fits to our JKL' 
photometry assuming pure hydrogen (top) and 
pure helium (bottom) atmospheres. The observed 
JK fluxes are represented by errorbars. Model 
monochromatic fluxes are shown as the solid red 
line. Our measurements yield physical parameter 
estimates as listed in Table 3. The L'-band non- 
detection from Keck is shown as an upper-limit. 



The transmission functions S'm(A) along with the 
constants Crr,. for each bandpass are de scribed at 



length in Holberg and Bergeronl ( 2006[ ) and ref- 
erences therein. To make use of all photomet- 
ric measurements simultaneously, we convert mea- 
sured magnitudes into observed fluxes using equa- 
tion (IlJ and compare the resulting energy distri- 
butions with those predicted from our model at- 
mosphere calculations. Model fluxes are also av- 
eraged over filter bandpasses by substituting f\ 
in equation ([2]) for the monochromatic Eddington 
flux H\. The average observed fluxes f\,m and 
model fluxes Hx^m, which depend on Tcff, log 5 
and iV(He)/iV(Hj, are related by 



/a,™ = 47T{R/dfH, 



(3) 



where R/d is the ratio of the radius of the star to 
its distance from Earth. 

Our fitting procedure relies on the nonlin- 
ear, least-squares, steepest descent method of 
Levenberg-Marquardt. The value of is taken 
as the sum over all bandpasses of the differ- 
ence between both sides of equation ([3]), prop- 
erly weighted by the corresponding observational 
uncertainties. In our fitting procedure, we con- 
sider only Toff and the solid angle, R/d, as free 



parameters. We begin by guessing a value for 
the surface gravity (logg = 8) and deter- 
mine the corresponding effective temperature and 
solid angle. Combined with the distance, c?, ob- 
tained from the trigonometric parallax measure- 
ment, this yields directly the radius of the star 
R. The radius is then converted into mass using 
evolutionary m odels similar to those described in 
Fontaine et al.l (jJoOl) but with C/0 cores, g(He) 
= log Mhc/M^ = 10-2 and q{H) = 10"^, which 
are representative of hydrogen atmosphere white 
dwarfs, and q{He) = 10"^ and qiH) = lO"!", 
which are representative of helium atmosphere 
white dwarfs. In general, the logg value obtained 
from the inferred mass and radius (g = GM/R'^) 
will be different from our initial guess of log g — 8, 
and the fitting procedure is thus repeated until an 
internal consistency in logt; is reached. 

Best fit results for the white dwarf physical 
properties are summarized in Table 3. Assum- 
ing a hydrogen-dominated atmosphere, we find: 
Tcff = 8160±4000 K and log 5 8.90±0.02, which 
corresponds to a cooling age of tc « 3.4 Gyr and 
model-dependent mass mmodci — 1-15 ± O.OIMq. 
Assuming a Helium-dominated atmosphere, we 
find: Tes = 7670 ± 4210 K, \ogg = 8.78 ± 0.02, 
tc ~ 3.5 Gyr, and Wmodoi = 1-08 ± O.OIMq. 
An m — 1.1 M0 remnant corres ponds to an 
6.2M(7) — 6.9M(n progeni tor mass fKalira i et al 



2008; W illiams et al]|2009l ). Thus, HD 114174 B 
spent only a short period of time (« 50 Myr) on 
the main-sequence. 

6. SUMMARY & DISCUSSION 

We report the discovery of a faint companion 
orbiting the nearby G5 star HD 114174. Initially 
identified as a promising high-contrast imaging 
target due to a measured long-term RV acceler- 
ation, we have used NIRC2 at Keck to obtain di- 
rect images of the source that is responsible for 
causing the Doppler trend. Multi-epoch detec- 
tions in two different filters demonstrate associ- 
ation through common proper- motion over an 1.4 
year time frame. The RV data spans 16.0 years 
and shows curvature (a change in the accelera- 
tion). Since HD 114174 has a precise parallax 
measurement, it will soon be possible to place 
strong constraints on the companion orbit (all six 
parameters) and mass as soon as the direct imag- 
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ing astrometry shows curvature. 

Our broadband photometric measurements in- 
dicate that the companion has neutral, or possi- 
bly blue colors {J — K = 0.12±0.16) suggesting 
interesting atmospheric physics. Deep ADI ob- 
servations in the L -band under good seeing con- 
ditions at Keck resulted in a non-detection, fur- 
ther supporting the interpretation of a decreasing 
NIR spectral-energy distribution with increasing 
wavelength. HD 114174 B has an absolute mag- 
nitude {Mj = 13.97 ±0.11) and colors consistent 
with a T3-dwarf. However, dynamical consider- 
ations preclude the interpretation of a substellar 
companion. We calculate a firm minimum mass 
of m > 0.260 ± O.OIOMq using presently available 
RV and imaging data. 

HD 114174 B appears to be a white dwarf. Us- 
ing theoretical models to fit our photometry and 
assuming a pure hydrogen atmosphere, we esti- 
mate Teff = 8160 ± 4000 K and \ogg = 8.90 ± 
0.02, which corresponds to a cooling age of tc « 
3.4 Gyr and remnant mass mmodoi — 1-15 ± 
O.OIM0, consistent with the 4.7l2'g Gyr age es- 
timate of the system based on isochronal anal- 
ysis of HD 114174 A, and also consistent with 
the dynamical mass lower limit. Calculations as- 
suming a Helium-rich atmosphere yield similar re- 
sults. Uncertainty in Toff is presently large be- 
cause the NIR colors of cool white dwarfs are 
known to exhibit an inversion (from having red 
colors to blue colors as T^s decreases) as the re- 
sult of collisionally-induced H2 absorption. Thus, 
there exists a degeneracy in Tes for J — if « 0.0. 
Deepe r L'-band observations obtained wi t h the 

Emit 



LBT (jSkrutskie et al 



Skemer et all |2012[ ) 



and direct NIR spect roscopy measuremen ts using 
an integral-field unit (jHinklev et al.ll201ir) will al- 
low us to refine HD 114174 B's physical proper- 
ties. Furthermore, it may be possible to study 
HD 114174 B at shorter wav elengths using HST in 
a fashion similar to that of Barstow et al.l ( 2005[ ) 
who obtained a spectrum of Sirius B with STIS. 

HD 114174 A is a marginally evolved Sun- 
like star located at d = 26.14 ± 0.37 pc. Thus, 
it may be possible to measure its radius di- 
rectly using CHARA interferometry, placing 
tighter constraints on the primary star's mass and 
hence the syste m age and also secondary mass 
( Bovaiian et al.l [2009). A similar analysis has 
been performed for the benchmark brown dwarf 



HR 7672 B (jCrepp et al.ll2012al) . In the case of a 
white dwarf, a more accurate mass and age will 
likewise constrain physical properties of the pro- 
genitor. The HD 114174 system thus represents a 
useful testbed for studying white dwarf initial-to- 
final mass ratios at a level of detail comparable to 
Sirius and Procyon. Given the anticipated sensi- 
tivity of forth-coming high-contrast imaging pro- 
grams that use "extreme" AO, it is likely that ad- 
ditional benchmark white dwarf companions will 
be (accidentally) uncovered in the near future. A 
more complete census of Sirius-like systems may 
ultimately be undertaken with the GAIA space 
mission using precision astrometry. 
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